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Abstract 

The oscillatory flow through tapered cylindrical tube sections (jet pumps) is characterized 
by a numerical parameter study. The shape of a jet pump results in asymmetric hydrody¬ 
namic end effects which cause a time-averaged pressure drop to occur under oscillatory flow 
conditions. Hence, jet pumps are used as streaming snppressors in closed-loop thermoacous¬ 
tic devices. A two-dimensional axisymmetric compntational flnid dynamics model is used 
to calculate the performance of a large nnmber of conical jet pump geometries in terms of 
time-averaged pressnre drop and acoustic power dissipation. The investigated geometrical 
parameters include the jet pump length, taper angle, waist diameter and waist curvatnre. In 
correspondence with previous work, four flow regimes are observed which characterize the jet 
pump performance and dimensionless parameters are introduced to scale the performance 
of the various jet pump geometries. The simulation results are compared to an existing 
qnasi-steady theory and it is shown that this theory is only applicable in a small operation 
region. Based on the scaling parameters, an optimum operation region is defined and design 
guidelines are proposed which can be directly used for future jet pump design. 


PACS numbers: 43.35.Ud, 43.25.Nm, 43.20.Mv, 47.32. Ff 
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FIG. 1: Jet pump with parameters that define the geometry (not to scale). Bottom dashed 
line indicates center line, top solid line indicates outer tube wall. Reproduced from^. 

I. INTRODUCTION 

One of the most promising configurations for thermoacoustic devices is the traveling wave 
configuration which generally consists of a closed-loop tube.^ This design however, has one 
disadvantage: the possibility of a time-averaged mass flow, or Gedeon streaming, to occur 
due to the looped geometrj^. This streaming leads to unwanted convective heat transport 
and reduces the efficiency of closed-loop thermoacoustic devices. 

A commonly used solution is the application of a jet pump.SEEEl This is a tapered 
tube section which, due to the asymmetry in the hydrodynamic end effects, establishes a 
time-averaged pressure drop. By balancing the time-averaged pressure drop across the jet 
pump with that which exists across the regenerator of a thermoacoustic device, the Gedeon 
streaming can be canceled.® Fig. [Tjshows a schematic of a typical conical jet pump geometry 
with its corresponding dimensions. The two openings both have a different radius: Rb for 
the big opening (left side) and Rg for the small opening. Together with the jet pump length, 
Ljp, the jet pump taper angle a is defined. Furthermore, at the small opening a curvature 
Rc is applied to increase the time-averaged pressure drop compared to a sharp contraction. 
To estimate the performance of a jet pump, a quasi-steady model has been proposed by 
Backhaus and Swift.® This model is based on minor losses in steady flow and assumes that 
the oscillatory flow can be approximated as two steady flows.® Given the pressure drop 
generated by a pipe transition in steady flow,® 

Apmi = ^Kpu‘^, ( 1 ) 

and assuming a pure sinusoidal velocity in the jet pump, the time-averaged pressure drop 
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across the jet pump can be estimated as 


^P2,JP 


gPo|Mi,jp| 



( 2 ) 


where is the velocity amplitude at the small exit of the jet pump. The subscripts “s” 

and “6” indicate the small and big opening of the jet pump, respectively. K^xp is the minor 
loss coefficient for abrupt expansion and i^con is the minor loss coefficient for contraction, 
both are well documented for steady flows.® 

With the jet pump being an additional flow resistance, acoustic power will be dissipated 
due to viscous dissipation and vortex formation. Using the same quasi-steady approach, 
Backhaus and Swift derived an equation to estimate the amount of acoustic power dissipation 
associated with the jet pump as® 


AUjp = 


Svr 


(^Kexp,s T Kcon,s) T 



2 

(^Kcon,b T Kexp,b^ 


(3) 


An optimal jet pump should establish the required amount of time-averaged pressure 
drop to cancel any Gedeon streaming with minimal acoustic power dissipation. This requires 
maximizing the difference in minor losses due to contraction and expansion while at the same 
time minimizing the sum of the minor loss coefficients. 

Previous studies have shown however, that the accuracy of this quasi-steady approx¬ 
imation is limited and that there are other factors influencing the performance of a jet 
pump.^ *^^ * ^^ * ^^ * ^^ ^ Petculescu and Wilen experimentally studied the influence of taper angle 
and curvature on minor loss coefficients.^ Good agreement was obtained between steady flow 
and oscillatory experiments for taper angles up to 10°. However, the minor loss coefficients 
determined were found to be strongly dependent on the taper angle used. A qualitative 


comparison between their hndings and the current results will be provided in Section III. B 


Smith and Swift investigated a single diameter transition, corresponding to one end of a 
jet pump.in^ The measured pressure drop and acoustic power dissipation was found to be 
dependent on the dimensionless stroke length, the dimensionless curvature and the acoustic 
Reynolds number. Recently, Tang et al. investigated the performance of jet pumps numer¬ 
ically,^ but assumed a priori the quasi-steady approximation to be valid by modeling the 
flow as two separate steady flows. The negative effect of flow separation on the jet pump 
performance was identihed which is in line with the current work. 

Although some jet pump measurements are available in literature, a systematic param¬ 
eter study that directly relates variations in wave amplitude and geometry to a jet pump’s 
performance has, to the authors knowledge, not yet been addressed. A hrst step towards 
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the investigation of a jet pnmp’s performance in oscillatory flows has been made in previous 
work by scaling the jet pump geometry using two different Keulegan-Carpenter numbers 
and correlating this to the jet pump performance.^ Four different flow regimes were distin¬ 
guished as a function of the wave amplitude. Fig. [^shows the observed vorticity helds in the 
vicinity of a jet pump at t = tmax when ujp{t) > 0. In order of increasing wave amplitude 
these flow regimes can be described as follows. At low wave amplitudes, an oscillatory vor¬ 
tex pair exists on both sides of the jet pump but no vortex shedding is observed (Fig. [^ . 
When the amplitude is increased, one-sided vortex propagation at the location of the small 
jet pump opening starts (Fig. 2b). In this flow regime, vortex rings are shed from the jet 
pump waist to the right while on the big side of the jet pump an oscillatory vortex pair is 
still visible. Vortex shedding from the big side of the jet pump begins to occur at increased 
wave amplitudes (Fig. 2c). Here vortex rings are shed from both jet pump openings directed 
outwards, but no interaction between the two sides is visible. Ultimately, a further increase 
in wave amplitude leads to interaction between the two jet pump openings (Fig. 2d). Vor¬ 
tices being shed from the small jet pump opening now propagate in alternating directions 
and flow separation inside the jet pump occurs. 

Albeit the described flow regimes were distinguished using the proposed scaling, the 
data set was too limited to determine whether the observed flow separation (Fig. 2d) was 
geometrically initiated by an increase in the taper angle or by a decrease in the jet pump 
length. In this paper, a parameter study is performed using a computational fluid dynamics 
(CFD) model to further identify the relation between the four different flow regimes and 
the jet pump geometry. The influence of various geometric parameters including the jet 
pump taper angle, length, curvature and waist diameter on the jet pump performance is 
investigated. Based on the presented results, the existing scaling parameters are further 
extended. Furthermore, a comparison with the quasi-steady model is provided to determine 
under what conditions the approximation is applicable. Based on this parameter study, 
design guidelines for future jet pump design are proposed. 

After a description of the used CFD model in Section [TT| the various investigated jet 
pump geometries are introduced in Section |II. D The resulting flow regimes will be dis¬ 
tinguished in Section |III. A| and subsequently linked to the jet pump performance in Sec¬ 
tion inri Finally, the influence of the jet pump curvature is investigated in Section nrn 


before drawing hnal conclusions on the design guidelines in Section IV 


II. MODELING 

A two-dimensional axisymmetric CFD model is developed using the commercial software 
package ANSYS CFX version 14.5,^ identical to the model used in previous work.^In the 
following, the numerical model with the corresponding boundary conditions will be repeated 
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(a) Oscillatory vortex pair on both sides and no jetting is observed. 



(b) Propagating vortex to right side and oscillating vortex pair on left side 
of the jet pump. 



(c) Propagating vortex on both sides of the jet pump. 



(d) Vortices propagate fully through jet pump and flow separation inside 
the jet pump occurs. 


FIG. 2: Four different flow regimes distinguished based on the instantaneous vorticity helds 
V X u [1/s] at t = tmax arouud the jet pump for the a = 7° geometry (color online). 
Reproduced froml^^. 
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briefly and the methods to derive the jet pump performance will be explained. 

The different jet pump geometries are conhned in an outer tube with radius Ro = 
30 mm. The length of the outer tube on either side of the jet pump is Lq = 500 mm when 
/ = 100 Hz. This length is scaled relative to the wavelength to avoid the jet pump being 
placed at a velocity node for simulations performed at other frequencies. In all cases, air at 
a mean pressure of po = 1 atm and a mean temperature of Tq = 300 K is used as the working 
fluid. 


A. Numerical model 


The unsteady, fully compressible Navier-Stokes equations are discretized using a high reso¬ 
lution advection scheme in space and a second order backward Euler scheme in time.^^ The 
system of equations is closed using the ideal gas law as the equation of state and the energy 
transport is calculated using the total energy equation, including viscous work terms. Based 
on the critical Reynolds number dehned in Section all presented results fall within the 
laminar regime so no additional turbulence modeling is applied. The time-step At is chosen 
such that each wave period is discretized using 1000 time-steps. A total of ten wave periods 
are simulated to achieve a steady periodic solution and the last hve wave periods are used 
for further analysis. 

The acoustic wave is generated on the left boundary of the domain using a sinusoidal 
velocity boundary condition with a specihed frequency and velocity amplitude. To control 
the wave propagation over the right boundary of the computational domain, a dedicated 
time-domain impedance boundary condition is developed and implemented in ANSYS CFX. 
The approach is based on the work of Polifke et ah and allows the application of a com¬ 
plex reflection coefficient at the boundary.^^^^ A detailed validation and explanation of the 
implementation can be found in the work of Van der Poel, which was carried out as part 
of the current research.^ In all cases described here, a reflection coefficient of |i?| = 0 is 
specihed to simulate a traveling wave on the right side of the jet pump. The combination of 
the velocity boundary condition and the time-domain impedance boundary condition results 
in a time-averaged volume how that is on average less than 0.5 % of the acoustic volume how 
rate. 


On the radial boundary of the outer tube (at r = i?o), a slip adiabatic wall boundary 
condition is imposed as the pipe losses in this part of the domain are currently not of interest. 
To correctly simulate the minor losses in the jet pump, a no-slip adiabatic wall boundary 
condition is used at the walls of the jet pump. 

The choice for a two-dimensional axisymmetric model to simulate how separation might 
need some additional explanation as in planar dihusers the how separation can be asymmet- 
ric | 20 | 21 | clear distinction should be made between how separation in planar and conical 
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diffusers. The steady flow separation in conical diffusers has been investigated previously 
and no visible asymmetry of the flow was reported under laminar conditions.^ Furthermore, 
the oscillatory nature of the flow prevents the flow from developing asymmetric instabili¬ 
ties, even in planar diffuser geometries.^^^^ These observations motivate the applicability 
of a two-dimensional axisymmetric model for the current situation of pure oscillatory flow 
through a conical geometry. 


B. Computational mesh 


The used spatial discretization is validated in previous work and will be only briefly described 
here.l^ The computational mesh consists of an unstructured part within 50 mm from the jet 
pump and a structured mesh in the rest of the domain. In both parts quadrilateral elements 
are used. In the jet pump region, a maximum element size of 1 mm is used which is rehned 
up to 0.5 mm near the jet pump waist. Moreover, to be able to accurately resolve the flow 
separation, a rehnement in the viscous boundary layer is applied such that a minimum of 
10 elements are located within one viscous penetration depth distance 6^ from the jet pump 
wall. For each simulated wave frequency, the mesh distribution is adjusted according to the 
viscous penetration depth, 5y = In the quadrilateral part of the mesh, a hxed 

radial mesh resolution of 1mm is used whereas the mesh size in the axial direction grows 
from 1 mm at a distance of 50 mm from the jet pump up to 10 mm at the extremities of the 
domain. 

In additional to the mesh validation performed in previous work,l^ the accuracy of 
the employed computational mesh has been investigated specihcally for the prediction of 
flow separation. The jet pump geometry with a taper angle of a = 15° (No. 4 in Table 0 
has been used at a frequency of 100 Hz and at a wave amplitude where flow separation is 
expected {^i/Dg ■ a = 0.85, Eq. [^. The mesh size in the jet pump region is in two steps 
rehned from 1 mm to 0.5 mm and 0.25 mm. At the same time, the mesh size in the jet pump 
waist region is rehned from 0.5 mm to 0.2 mm and 0.1 mm. Additionally, the ehect of a mesh 
rehnement in the viscous boundary layer is investigated by increasing the number of elements 


from 10 to 20. The diherence in the dimensionless pressure drop (Eq. 10) and dimensionless 


acoustic power dissipation (Eq. 11) between the investigated meshes is less than 0.04 and 
0.05, respectively. Furthermore, both the time and location where the how hrst separates 
do not vary more than At/T = 8- 10“^ and Ax = 0.4 mm between the diherent mesh 
resolutions. 
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C. Data analysis 

In order to determine the jet pnmp performance from the transient CFD solution, some ad¬ 
ditional analysis steps are performed. The first order amplitudes of all physical quantities are 
calculated by a discrete Fourier transformation for the specihed wave frequency using the so¬ 
lution from the last five simulated wave periods. To obtain the time-averaged quantities, the 
CFD solution is averaged over the last five wave periods in order to eliminate the first order 
contribution. The streaming velocity is calculated using a density-weighted time-average: 
U 2 = (pu)/(p). The velocity amplitude in the jet pump waist, |mi,jp|, is calculated using an 
area-weighted average over the cross-section. Subsequently, the local acoustic displacement 
amplitude in the jet pump waist can be calculated under the assumption of a sinusoidal 
jet pump velocity: 

^ 27r/ 

To scale the displacement amplitude to the jet pump geometry, a Keulegan-Carpenter num¬ 
ber is defined based on the jet pump waist diameter: 

KCd = F. (5) 

^ S 

where Dg = 2Rs is the jet pump waist diameter. This KCd is also referred to as a “dimen¬ 
sionless stroke length” or “dimensionless displacement amplitude” and can be rewritten to 
Rcd/S^ where Rcd is an acoustic Reynolds number based on diameter and S = 
the Stokes numberl^ ^ * ^^ * ^^ l 

D. Jet pump geometries 

The described computational model is used to investigate a large variety of conical jet pump 
geometries and wave amplitudes. By varying either the jet pump length, the size of one of 
the two openings or the local curvature at the jet pump waist, 19 different geometries are 
considered. The different investigated geometries are all variations based on one reference 
geometry (denoted by “re/” in all tables). The reference geometry is a jet pump with a 7° 
taper angle, all other reference dimensions are denoted in Table 

Varying the jet pump taper angle while keeping the waist diameter constant can be 
achieved by either changing the jet pump length Ljp or changing the radius of the big 
opening Rb. Because the jet pump length is not taken into account in the quasi-steady 
approximation (Eq. |^, the predicted jet pump performance is independent of the jet pump 
length. In a similar way, the radius of the big opening is only a weak parameter in the 
quasi-steady approximation as long as Rb 3> Rg- Both approaches are considered here 
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TABLE I; Dimensions of the reference jet pump geometry. 

Rq 30 mm 

Rh 15 mm 

Rs 7 mm 

Rc 5 mm 

Ljp 70.5 mm 
a 7“ 

TABLE If; Dimensions of jet pump geometries with varied taper angles a by either changing 
the jet pump length or the big radius. Constant radius of the jet pump waist Rg = 7.0 mm 
and constant curvature of the small opening Rc = 5.0 mm. 


No. 

symbol 

a 

Ljp 

Rb 

1 

• 

3° 

157.8 mm 

15.0 mm 

2 

■ 

5° 

96.7 mm 

15.0 mm 

ref 

X 

70 

70.5 mm 

15.0 mm 

3 

T 

0 

0 

50.8 mm 

15.0 mm 

4 

♦ 

15° 

35.5 mm 

15.0 mm 

5 

► 

20° 

27.9 mm 

15.0 mm 

6 

• 

3° 

70.5 mm 

10.4mm 

7 

■ 

5° 

70.5 mm 

12.7 mm 

8 

T 

0 

0 

70.5 mm 

18.5 mm 

9 

♦ 

15° 

70.5 mm 

24.4 mm 


to see what actually influences the jet pump performance: the distance between the two 
openings (Ljp) or the taper angle of the jet pump inner surface (a). Table [ll| shows the 
dimensions of the investigated geometries where the taper angle is changed with respect to 
the reference case by either varying the jet pump length (cases 1-5) or by varying the big 
radius (cases 6-9). The jet pump waist radius and curvature are kept constant at = 7 mm 
and Rc = 5 mm, respectively. This results in a dimensionless curvature of y = Rc/Dg = 0.36 
which is considered a “smooth” contraction.® Correspondingly, the contraction minor loss 
coefficient at the jet pump waist {Kcon,s in Eq. is expected to be negligible which will 
enhance the overall time-averaged pressure drop. 

In a different set of simulations, both the taper angle and the jet pump length are hxed 
while the jet pump waist diameter is changed (cases 10-13). Furthermore, the jet pump 
curvature is adjusted to maintain the same dimensionless curvature. The dimensions of this 


set of geometries are shown in Table III 
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TABLE III; Dimensions of jet pump geometries with varied waist radius Rg. Constant taper 
angle a = 7° and jet pump length Ljp = 70.5 mm. 


No. 

Rs 

Rb 

Rc 

10 

3.0 mm 

11.4 mm 

2.1 mm 

11 

5.0 mm 

13.2 mm 

3.6 mm 

ref 

7.0 mm 

15.0 mm 

5.0 mm 

12 

9.0 mm 

16.8 mm 

6.4 mm 

13 

11.0 mm 

18.6 mm 

7.9 mm 


In the following, the observed flow regimes and jet pump performance will be discussed 
for the cases where either the jet pump length or the size of one of the two openings is varied. 
In Section [III. C[ the effect of the jet pump waist curvature will be investigated. 


III. RESULTS AND DISCUSSION 

For each geometry, various wave amplitudes and consequently, various displacement ampli¬ 
tudes are simulated. The wave amplitude is varied by specifying different values of the ve¬ 
locity amplitude at the left boundary condition. It is important to ensure that the simulated 
wave amplitudes fall within the laminar regime as URANS turbulence models, especially 
in combination with wall functions, are known to lose their validity for relaminarizing and 
separating flows which is an important part of the current data set.l ^°*^^ * ^ ^ Other approaches 
to model turbulence such as large eddy simulation (LES) or direct numerical simulation 
(DNS) are currently not feasible for a parameter study of the size presented here. In order 
to determine what wave amplitudes fall within the laminar regime, the critical Reynolds 
number derived by Ohmi and Iguchi for oscillatory pipe flow is used,*^^*^ 

Re, = 305 ' , (6) 

where the Reynolds number is defined based on the viscous penetration depth; Re = 
\ui\5ypo/ hq. Normally, the velocity amplitude is directly proportional to the frequency for a 
given displacement amplitude, but the velocity amplitude where a transition to turbulence 
occurs {Re = Re,) scales with From previous work, it was concluded that the jet pump 
performance scales with the acoustic displacement amplitude instead of the jet pump velocity 
amplitude.^ Consequently, by lowering the driving frequency all displacement amplitudes 
of interest are investigated. The majority of the presented simulations are carried out with 
/ = 100 Hz while some additional simulations are using frequencies of 10 Hz, 20 Hz, 30 Hz, 
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50 Hz and 200 Hz. All data presented in this paper represents a total of 197 simnlations 
having a total single core compntational time of 391 days on an Intel Core i7 CPU. 

From the simnlated jet pnmp performance, a scaling parameter is introdnced which 
best aligns the jet pnmp performance for the varions investigated geometries and correctly 
incorporates the effect of the jet pnmp taper angle: 



(7) 


This scaling parameter is essentially the Kenlegan-Carpenter nnmber (Eq. mnltiplied by 
the taper angle (in radians). 


A. Flow regimes 


Becanse the jet pnmp performance is fonnd to be a strong fnnction of the flow regime, the 
observed flow regimes will hrst be discnssed. By stndying the transient vorticity held and 
streaming velocity held of all simnlations, a diherentiation is made between the fonr how 
regimes that have been illustrated in Fig. The how regimes that have been observed using 
various wave amplitudes and jet pump geometries are all presented in the {^i/Dg ■ a, KCd) 
space in Fig. 3a As such, the horizontal axis includes the taper angle while the vertical 


axis represents the Kenlegan-Carpenter number without the contribution of the taper angle. 
The diherent symbols represent the diherent how regimes: circles for oscillatory vortex pairs, 
squares for right-sided vortex propagation, triangles for two-sided vortex propagation and 
diamonds for how separation and alternating vortex shedding. The open circles indicate the 
position where the maximum dimensionless pressure drop was achieved with a specihc jet 
pump geometry. 

Two how regimes are easily distinguished in Fig. For small values of the Kenlegan- 
Carpenter number, KCu < 0.5, no vortex shedding is observed. This is indicated by the 
horizontal dashed line in Fig. In this regime, no asymmetry in minor losses occurs and 
consequently the time-averaged pressure drop is zero. The transition from this how regime 
to the situation where vortex shedding is hrst observed is independent of the taper angle 
and occurs when KCd > 0.5. Above this value the how held becomes asymmetric resulting 
in a positive time-averaged pressure drop. This is in line with the work of Holman et ah 
on synthetic jets where a formation criterion for vortex shedding is derived of the form 
KCd = Rcd/S'^ > C where C is a constant depending on the geometry of the orihce.^ 

In Fig. I^the how separation regime (diamonds) can be distinguished from the other 
three how regimes and occurs when ^i/Dg ■ a > 0.7 (dashed vertical line). Because this 
transition is clearly determined by the scaling parameter as introduced in Eq. the taper 
angle plays an important role in the initiation of the how separation regime. This how regime 
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is dominated by vortices being shed from the jet pump waist in alternating directions. This 
results in the flow held again becoming more and more symmetric. Despite the interesting 
how phenomena occurring in this how regime,^ it is of little practical interest when designing 
jet pumps as the achieved jet pump ehectiveness is low. Consequently, the area of interest 
for jet pump applications is bounded to KCd > 0.5 and ii/Dg ■ a < 0.7. 

Inside this area, two how regimes occur which cannot solely be separated based on KCd 
or ^i/Dg- a. Whether vortex shedding from the big jet pump opening occurs, is found to be 
dependent on the the local displacement amplitude at the position of the big opening 
A Keulegan-Carpenter number at the big jet pump opening is estimated using 


C -^s 

KCnfi = 

Uh 


K 

2Rl 


( 8 ) 


Note that is still the local displacement amplitude at the jet pump waist. Fig. [3b] shows 
the how regimes where either one-sided or two-sided vortex propagation occurs with KCofi 
instead oi KCd on the vertical axis. Only the range 0 < ^i/Dg-a < 0.7 is shown, omitting the 
how separation regime. The two-sided vortex propagation is hrst observed when KCD,b > 
0.15. This is at a diherent value than where vortex propagation at the jet pump waist 
was hrst recognized {KCd > 0.5), which is an ehect of the local dimensionless curvature 
y. The observed transition at KCD,b = 0.15 very nicely corresponds to the value found 
experimentally for a sharp-edged axisymmetric orihce.^^The ehect of the jet pump curvature 


will be further investigated in Section III. C 


B. Jet pump performance 

Ideally, a jet pump should impose the desired time-averaged pressure drop while having 
minimum acoustic power dissipation. Thus, the jet pump performance can be determined 
using two quantities: the time-averaged pressure drop Ap 2 and the acoustic power dissipation 
Ai725 both a function of the jet pump waist velocity amplitude |mi,jp|. The time-averaged 
pressure drop over the jet pump is calculated from the simulation results by a subtraction 
of the mean pressure helds on either side of the jet pump. In this way the inhuence of the 
vortex propagation distance on the time-averaged pressure drop is eliminated.!^ 

The acoustic power is calculated using hrst order quantities as®! 

E2{x) = ]^iR\pi{x)Ui{x)], (9) 

with pi{x) the complex conjugate of the pressure amplitude area-averaged over the local 
cross section and Ui{x) the complex volume how rate. The dissipation of acoustic power 
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(a) All four flow regimes. Dashed lines determine the(b) Flow regimes with one-sided (■) and two-sided 
bounds of the different flow regimes. (►) vortex propagation. The other two flow regimes 

are not shown. The local Keulegan-Carpenter number 
at the big jet pump opening, KCo,b, is on the vertical 
axis. The flow regimes are separated by horizontal 
dashed line at KCD,b = 0.15. 


FIG. 3: Observed flow regimes using jet pump geometries with varied taper angle and waist 
diameter plotted in two variable spaces. The different flow regimes are distinguished using 
different symbols: • for oscillatory vortices, ■ for right-side vortex propagation only, ► 
for two-sided vortex propagation and ♦ for flow separation. The dashed lines represent the 
determined bounds of the flow regimes and the open circles indicate the points where the 
maximum Ap 2 was found for a specific jet pump geometry (color online). 
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across the jet pump, AE 2 , is determined in a similar manner as the time-averaged pressure 
drop. 

To better understand the relation between these two quantities, a scaling is applied 
according to the work of Smith and Swift 


ae; = 


8Ap2 

Pol^i, jpP ’ 
811AE2 

pQ'nRl\uipp 


3 ■ 


( 10 ) 

( 11 ) 


Note that under the quasi-steady approximation, Ap^ would represent the difference in minor 
loss coefficients between the two ffow directions while AE 2 would represent the summation 
of the minor loss coefficients (i.e. the last terms in Eq. [^and[^ respectively). Using these 
two dimensionless quantities, a jet pump effectiveness can be defined;!^ 


Ap^ 

AEl 


( 12 ) 


Optimal jet pump performance is achieved when p is maximized. These three performance 
characteristics are calculated for each simulation case. In the following, they will be shown 
as a function of the scaling parameters introduced previously to study the effect of geometry 
changes on the jet pump performance. 


1. Time-averaged pressure drop 


Fig. 4a shows the dimensionless pressure drop Ap^ as a function ii/Dg ■ a. The symbols 


represent different taper angles, obtained by either changing the jet pump length or the big 
radius, each depicted in Table The results with a varied waist diameter are all plotted 
using crosses (x), representing their taper angle of a = 7°. The maximum dimensionless 
pressure drop for all cases is achieved between 0.2 < ^i/Dg-a < 0.5 which fits well within the 
determined range of the flow regimes where vortex shedding, but no flow separation, occurs. 
This emphasizes the direct relation between the flow field asymmetry and the time-averaged 
pressure drop. 

A decay in Ap 2 is observed for higher values of ^i/Dg ■ a. The range where flow 
separation was observed coincides with this decay in dimensionless pressure drop. The decay 
follows the same trend for all jet pump geometries and is properly scaled using the adjusted 
Keulegan-Carpenter number, which shows the dependency of the flow separation occurrence 
on the taper angle. Moreover, a corresponding decay in Ap 2 as a function of ^i/Dg ■ a is 
observed in preliminary experimental results that have been obtained at Re < Rcc using the 
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(a) Dimensionless pressure drop as a function of ^i/Ds-(b) Dimensionless pressure drop as a function oi KCd, 
a only the range KCd < 3 is shown. 

FIG. 4: Dimensionless pressure drop using jet pump geometries with varied taper angle 
(Table 0 or waist diameter (Table EH- The symbols represent different taper angles as 
listed in Table The thin dashed lines indicate the upper and lower values of Ap 2 from the 
quasi-steady approximation (Eq. 1^ for the various geometries (color online). 


15° jet pump geometry in an oscillatory flow facility.l^IESI By performing hot-wire anemometry 
near the big jet pump opening the existence of an outward directed vortex street was revealed 
when ^i/Dg ■ a > 0.7, confirming the observations from the current simulations. A detailed 
discussion of the experimental results is deferred to a future publication. 

In the published data of Petculescu and Wilen however, the decay in dimensionless 
pressure drop is less severe.^ Their values of Ap 2 typically stabilize between 0.20 and 0.50 
at high values of ^i/Dg ■ a, depending on the taper angle used, and do not decrease to zero. 
Two possible explanations will be discussed: their small waist diameter compared to the 
viscous penetration depth and a transition to turbulence. In the geometries of Petculescu 
and Wilen, ~ 8.3 while for the current geometries ranges from 27 to 98. When 

the waist diameter is on the order of the viscous penetration depth, the occurrence of flow 
separation, and consequently the dimensionless pressure drop at ^i/Dg-a > 0.7, is most likely 
affected due to viscous effects dominating the flow and preventing the flow from separating. 
Alternatively, turbulence can affect the decay in Ap 2 when the turbulent boundary layer 
prevents the flow from separating.®' A detailed investigation of the effect of turbulence on 
the oscillatory flow in jet pumps is subject to further research. 

The adjusted Keulegan-Carpenter number, ^i/Dg- a, provides consistent scaling of the 
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dimensionless pressure drop in the decaying part of the curve shown in Fig. 4a However, a 
large spread from the reference geometry is observed for low values of the adjusted Keulegan- 
Carpenter number {^i/Dg ■ a < 0.2). Fig. 4b shows the dimensionless pressure drop as a 
function of KCe> which better scales the increase of Ap 2 towards its maximum. This shows 
that the initial increase of Ap 2 towards its maximum is not dependent on the taper angle 
at all. A clear similarity in this graph is the range KCu < 0.5 where the pressure drop is 
negligible, corresponding to the flow regime with oscillatory vortices shown in Fig. At 
larger values of KCd, the effect of the taper angle becomes apparent and the results of the 
different geometries deviate from each other. 


2. Acoustic power dissipation 

By studying the dimensionless acoustic power dissipation, it is observed that the jet pump 
taper angle has less effect on this performance quantity compared to the pressure drop. 
Instead, the dimensionless acoustic power dissipation (Ai?^) shows good similarity between 
the different jet pump geometries when scaled using KCd- Fig. shows the dimensionless 
acoustic power dissipation as a function of the Keulegan-Carpenter number based on the jet 
pump waist diameter, KCd- The symbols correspond to those in Fig.|^ A common minimum 
occurs around KCd ~ 1 while at higher amplitudes {KCd > 3) AE 2 only slightly increases. 
Although the taper angle has no effect on the location of the minimum dimensionless acoustic 
power dissipation, the slight increase of AE 2 at higher amplitudes is easily observed for higher 
taper angles (15°, ♦ and 20°, ►). 

3. Comparison to quasi-steady model 

The obtained simulation results can all be compared to the quasi-steady model proposed 
by Backhaus and Swift. SI In this model, the values of Ap 2 and AE^ are assumed to be 
amplitude independent and the only signihcant variation in the quasi-steady model occurs 
by changing the waist diameter. In Fig. |4]-[^ the upper and lower bounds of these quasi¬ 
steady approximated values are indicated by the thin dashed lines. It becomes clear that the 
quasi-steady approximation is a good approach to estimate the time-averaged pressure drop 
of a jet pump when operating in the optimum ranges of ii/Dg - a and KCd- Outside of these 
ranges however, the actual performance will deviate signihcantly from the predicted values. 
The simulated acoustic power dissipation is especially in good agreement for the lower taper 
angles and KCd > 2. For high taper angles (15°, ♦ and 20°, ►) however, a deviation from 
the quasi-steady model is observed at KCd > 2. On average the quasi-steady values show 
considerable correspondence to the measured dimensionless acoustic power dissipation. 
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FIG. 5: Dimensionless acoustic power dissipation using jet pump geometries with varied 
taper angle (Table ig or waist diameter (Table |nT|). Symbols in accordance with Fig. || 
(color online). The thin dashed lines indicate upper and lower values of the predicted acoustic 
power dissipation according to the quasi-steady approximation (Eq. [^. 


4- Jet pump effectiveness 

The two jet pump performance characteristics, dimensionless pressure drop and acoustic 
power dissipation, can be combined into one effectiveness parameter p as dehned in Eq. [T^ 
Fig. 1^ shows the jet pump effectiveness for all geometries with varied taper angle or waist 
diameter. The highest jet pump effectiveness is again obtained between 0.2 < f^i/Dg-a < 0.5. 
Note that this does not strictly coincide with the minimum in the acoustic power dissipation 
(at KCd ~ !)• However, as long as KCd > 0.5, the variations in dimensionless pressure 
drop have a larger influence on the jet pump effectiveness than the relative small changes in 
the dimensionless acoustic power dissipation. 


Summarizing the results so far, it can be concluded that for laminar oscillatory air 
flows the optimum jet pump performance is achieved between 0.2 < f,i/Dg ■ a < 0.5 under 
the condition that KCd > 0.5. The combination of these two criteria imposes a limit on 
the maximum taper angle which maximizes the dimensionless pressure drop. Substitution 
of the two criteria, yields a < 23° in order to prevent overlapping of the range of optimum 
jet pump performance with the range where no pressure drop is measured. The optimum 
jet pump performance criterion does well coincide with the region where either right-sided 
or two-sided vortex shedding is observed (indicated with ■ and ► in Fig. 3a). It is only in 
this optimum region that the measured dimensionless pressure drop will approach values as 
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FIG. 6: Jet pump effectiveness for geometries with varied taper angle (Table or waist 
diameter (Table [T^. Symbols in accordance with Fig. (color online). 


predicted by the quasi-steady model. 

C. Influence of waist curvature 

In the previous discussion where the influence of the jet pump taper angle and waist diameter 
were investigated, the radius of curvature at the jet pump waist was kept constant. In steady 
flow theory, the minor loss coefficient for contraction is dependent on the dimensionless 
radius of curvature,® y = Rc/Dg. In order to maximize the time-averaged pressure drop 
while minimizing the acoustic power dissipation, it is desirable to minimize the contraction 
minor loss coefficient (following Eq. |2]-[^ and as such, apply a smooth curvature at the jet 
pump waist. For a “smooth” contraction (y > 0.15) the theoretical minor loss coefficient is 
Kcon = 0.04 which increases up to Kcon = 0.50 for a sharp contraction.® For all the jet pump 
geometries shown so far, a curvature of y = 0.36 was applied. In the current section, the jet 
pump curvature will be varied to study both its effect on the jet pump performance as well 
as on the occurring flow phenomena. 

Six different dimensionless curvatures have been investigated, ranging from y = 0 to 
0.65. Aside from the radius of curvature i?c, the jet pump geometry is equal to the reference 
geometry found in Table [T| An overview of the studied geometries with variable curvature is 
shown in Table HYl 
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(c) 

FIG. 7: Dimensionless pressure drop using jet pump geometries with varied waist curvature 
as a function of three different scaling parameters: a) KCc, b) iFC'£)/0.15(l + x)‘^ and c) 
KCd ■ (1 + 0.15/x)- The different symbols correspond to dimensionless curvatures ranging 
from y = 0 to 0.65. The dashed black lines indicate the upper and lower values from the 
quasi-steady approximation for the various values of y (color online). 
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TABLE IV: Dimensions of jet pump geometries with varied radius of curvature Rc. Constant 
taper angle a = 7°. Horizontal line separates cases at x = 0.15. 


No. 

Rc 

X 


14 

0.0 mm 

0 

“sharp” 

15 

1.0 mm 

0.08 

16 

2.5 mm 

0.19 


ref 

5.0 mm 

0.36 

“smooth” 

17 

7.0 mm 

0.48 

18 

10.0 mm 

0.65 



1. Jet pump performance 


When comparing the measured dimensionless pressure drop using these jet pump samples, 
it is found that an additional contribution to the Keulegan-Carpenter number is required to 


properly scale the results from different curvatures. Fig. [7a] shows the dimensionless pressure 
drop as a function of KCd for the geometries with varied curvature. The investigated 
curvatures are indicated by different symbols depicted in the corresponding legend. The 
crosses (x) represent the reference geometry where y = 0.36. The shape of the dimensionless 
pressure drop curve is found to be inversely proportional to the dimensionless curvature, 
especially for low values of y. For large curvatures, the effect on the dimensionless pressure 
drop becomes less apparent. Holman et al. studied the effect of edge curvature on the onset 
of a synthetic jet.^^ A formation criterion of the form KC^ > C ■ {1 + xY was derived 
based on the vorticity flux through the jet. The applicability of this criterion to jet pumps 
with variable curvature is investigated in the following. It is observed that especially the 
initial increase in Ap 2 , corresponding to the onset of vortex shedding, can be properly scaled 
using KCd/C{1 + xY with p = 4. Given that the onset of vortex shedding for a jet pump 
geometry with y = 0.36 is found at KCd = 0.5, the additional scaling parameter C = 0.15 
can be estimated. Fig. [7b] shows the dimensionless pressure drop as a function of the adjusted 
scaling parameter. While the initial increase in Ap 2 for the various curvatures all collapse 
to a single curve, the decay in Ap 2 due to the occurrence of flow separation cannot properly 
be accounted for by any scaling of the form KCd/C{1 + xY- 

To overcome this, an alternative scaling strategy for y is proposed of the form KCd • 
(1 + a/y). Using a least-squares approach, the htting parameter a = 0.15 is determined. 
Note that 0.15/y represents the fraction of the dimensionless curvature with respect to a 
“smooth” contraction. Fig. [7c| shows again the dimensionless pressure drop data, but now 
plotted against KCd ■ (1 + 0.15/y). The case with no curvature (y = 0) does not appear in 
this hgure as it will lead to inhnite values on the horizontal axis. Both the maxima of the 
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dimensionless pressure drop as well as the decay at higher amplitudes are well aligned using 
this alternative scaling approach. 

What is remarkable is that the maximum measured dimensionless pressure drop de¬ 
creases with increasing curvature, which is in contrast with the steady flow theory where 
the minor loss coefficient for contraction is expected to fall until y > 0.15. The maximum 
achieved dimensionless pressure drop ranges from Ap^ = 1.20 for a sharp contraction (y = 0) 
to Ap 2 = 0.69 for y = 0.65 in an almost linear manner. Furthermore, the effect of curvature 
on the dimensionless pressure drop is not only reversed, but the variation is larger than what 
is predicted by the quasi-steady model (dashed lines). An explanation of this behavior is 
the effect of curvature on the expansion phase, which is typically neglected in a quasi-steady 
approximation of the jet pump performance. Not only will the minor loss coefficient for 
contraction at the jet pump waist decrease, but so does the minor loss coefficient for expan¬ 
sion. This was first experimentally determined by Smith under steady flow conditions.®' A 
decrease in the expansion minor loss coefficient at the jet pump waist will lead to a decrease 
in the dimensionless pressure drop according to Eq. 

Additionally, this is reflected in the measured dimensionless acoustic power dissipation. 


which is shown in Fig. ^ The minimum achieved dimensionless acoustic power dissipation 
decreases more than six times between y = 0 (not shown) and y = 0.65. A similar effect 
has been reported when rounding the orifice of an synthetic jet.®' Combining the effect 
the curvature has on both the pressure drop and acoustic power dissipation using the jet 
pump effectiveness, shows that a higher dimensionless curvature leads to a better jet pump 
performance. This is depicted in Fig. |8b| where the jet pump effectiveness for five investigated 
curvatures is shown. 


2. Flow regimes 


Similar to what was presented in Section III. A the observed jet pump behavior can be ex¬ 
plained by studying the flow field and distinguishing between the four different flow regimes. 
Fig-H shows the distribution of the flow regimes in the (y, KCd) space. The four flow 
regimes are represented by different symbols in accordance to the symbols used in Fig. A 
clear influence of the dimensionless curvature is visible on the value oi KCd where vortex 
shedding is first observed (transition from • to ■) as well as on the value where flow sepa¬ 
ration is initiated (transition from ► to ♦). The boundary between one-sided and two-sided 
vortex shedding (transition from ■ to ►) is not determined by the jet pump waist curvature 
at all. As the local curvature at the big opening is not changed, no dependency on y is to 
be expected. 

By substituting the original defined flow regime bounds from Section III. A| in the ad¬ 
justed scaling parameter KCo ■ (1 -|-0.15/x), the bounds between the flow regimes (indicated 
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FIG. 8: Dimensionless acoustic power dissipation (a) and jet pump effectiveness (b) using 
jet pump geometries with varied waist curvature and scaled using KCo ■ (1 + 0.15/x)- The 
different symbols correspond to dimensionless curvatures ranging from y = 0.08 to 0.65. The 
dashed black lines indicate the upper and lower values from the quasi-steady approximation 
for the various values of y (color online). 
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by symbols as used in Fig. |^and Fig. can be represented in the (y, KCn) space: 


• ^ KCd = 0.7 (^1 + 

d3 

■ ^ KCd = 0.15^, 


0.15 

X 


-1 


A 1 0-15 

—)■ ♦: KCd — “ ( 1 “I- 

« V X 


-1 


(13) 

(14) 

(15) 


where the hrst expression corresponds to KCd = 0.5, the second expression corresponds 
to KCD,b = 0.15 and the last expression corresponds to ii/Dg ■ a = 0.7 for the cases 
where y = 0.36. The bounds are shown in Fig. ^ by the dashed lines. The bottom curve 
represents the initiation of vortex shedding (Eq. 13) and the upper curve the transition to 


flow separation (Eq. 15), both a function of y. The horizontal dashed line in Fig. [^indicates 
the boundary between one-sided and two-sided vortex shedding. 

Based on the data available the determined bounds do separate the flow regimes de¬ 
picted in Fig. I^well, although the defined onset of vortex shedding (• —■) is subject to 
discussion. Given the fact that the initial increase in Ap 2 caused by the onset of vortex 
shedding was well described by adopting a scaling of the form proposed by Holman et al.,1^ 
a flow regime bound based on the latter could be more appropriate: 


• ^ 


A'Cd = 0.15(1+ ,\f 


(16) 


This is represented in Fig. by a gray solid line and does indeed separate the two flow 
regimes. It should be noted that the data available is not sufficient to determine whether 


Eq. 13 or Eq. 16 is more appropriate. 


Nevertheless, Eq. IT predicts vortex shedding for y; = 0 at KCd > 0.15 which is 
also the formation criterion observed experimentally for a sharp edged synthetic jet.^ This 
exactly describes the difference in the transitional Keulegan-Carpenter number, observed 


in Section UFA, between vortex shedding from the jet pump waist {KCd,s > 0.5) and 
vortex shedding from the big jet pump opening {KCD,b > 0.15). A lower dimensionless 
curvature leads to a lower value of KCd where vortex shedding is initiated. Consequently, 
vortex shedding at the sharp edged big opening will take place at a lower value of the local 
Keulegan-Carpenter number than the vortex shedding from the jet pump waist. 


IV. CONCLUSIONS 

A CFD based parameter study is performed to investigate the influence of various geometric 
parameters on the performance of jet pumps for thermoacoustic applications. A total of 19 
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FIG. 9: Distribution of the different flow regimes in the {x, KCd) space for jet pump 
geometries with varied curvature (Table [Iv]) and a 7° taper angle. Symbols are in accordance 
with Fig. with the open circles showing the position where Ap 2 is maximum for each 
investigated curvature. The dashed lines represent the bounds of the different flow regimes 
determined according to Eq. lA -lA and the gray solid line represents the onset of vortex 
shedding according to Eq. 16 (color online). 


different jet pump geometries are used and for each geometry a number of wave amplitudes 
are simulated, resulting in a total of 197 simulations. In correspondence with previous 
work, four flow regimes are distinguished and separated in a hxed variable space.^ This 
space spans the jet pump taper angle a, the dimensionless curvature y and the Keulegan- 
Carpenter number KCd. At a certain value of KCd, single-sided vortex propagation from 
the jet pump waist is initiated and the flow held becomes asymmetric between the left and 
right side of the jet pump. This onset is found at KCd = 0.5 for a dimensionless curvature 
of y = 0.36. Two different expressions are proposed to account for the effect of y, but more 
data is required to decisively choose one of the two. At some point, how separation inside the 
jet pump can be distinguished and vortices are shed from the jet pump waist in alternating 
directions leading to a how held that is more symmetric. This transition is found to be 
strongly dependent on the jet pump taper angle and occurs when KCd > ^(l + 0-15/y)~\ 
The amount of asymmetry in the how held has a direct consequence on the mea¬ 
sured time-averaged pressure drop. When no vortex shedding occurs, no dimensionless time- 
averaged pressure drop is observed. Furthermore, as soon as the how separation is observed, 
the dimensionless pressure drop has already decayed signihcantly compared to the maximum 
value. For all the performed simulations, the maximum jet pump ehectiveness is observed 










Oosterhuis, Biihler, Wilcox and Van der Meer, JASA, p. 26 


when either single-sided or two-sided vortex propagation occurs. Consequently, the practical 
operation area of jet pumps with Dg/Sy 3> 1 and Re < Rcc is bound between the onset of 
vortex shedding from the jet pump waist and the occurrence of flow separation. This can be 
used as a guideline for future jet pump design. 

The design considerations presented in this paper, are supposed to provide better insight 
into the validity of the quasi-steady approximation that is widely used for the design of jet 
pumps for thermoacoustic applications. The quasi-steady approximation is in most cases 
an ideal representation of the jet pump performance and is valid only in a small operation 
area. Outside this region, the actual jet pump performance is expected to be signihcantly 
lower and a correction using the presented results is advised. 

Although the current parameter study represents a wide range of jet pump geometries, 
it is restricted to jet pumps with a single, linear tapered hole. As such, this work can serve 
as a basis for further jet pump geometry optimization in terms of effectiveness, compactness 
and robustness. Extensions to multiple holes or different shaped jet pump walls might affect 
the jet pump performance and the occurrence of flow separation. Furthermore, the influence 
of turbulence on the flow separation requires further attention. If the flow separation regime 
could be shifted to higher Keulegan-Carpenter numbers without a decrease in effectiveness, 
the operation area of a jet pump can be extended, greatly enhancing its robustness. 
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